SECURITY  CLASSIFICATION  OF  THIS  FAGS 


CONTENTS 


INTRODUCTION  .  1 

THEORY  OF  CONVECTIVE  COOLING  .  2 

NUMERICAL  SIMULATIONS  .  6 

CONCLUSIONS .  9 

ACKNOWLEDGMENTS .  9 

REFERENCES  . 14 


DTIC 

SELECTE|% 

APR2  W«JJ 


Accession  For 

/ 

NTIS  GRA&I 

DTIC  TAB 

□ 

Unannounced 

□ 

Justification _ 

Bv 

Distribution/ 

Availability  Codes 


THEORY  OF  BEAM  CHANNEL  HYDRODYNAMICS 


Introduction 


The  propagation  of  an  intense  beam  through  a  background  gas  can  result 
in  the  deposition  of  significant  energy  along  the  beam  path,  producing  a  hot, 
reduced-density  channel,  the  properties  of  the  channel  will,  in  turn,  affect 
subsequent  beam  pulses  which  travel  along  the  same  path.  In  recent 
experiments,  Greig  et  al.  [1-3]  have  studied  the  evolution  of  channels 
produced  by  lasers  and  electric  discharges  in  air  or  nitrogen  at  approximately 
1  atm.  From  those  data  a  general  picture  of  the  dynamics  of  beam  channels  has 
emerged,  including  time  scales  characterizing  channel  evolution. 

In  this  paper,  we  will  consider  pulses  which  deposit  energy  on  time  scales 
which  are  short  compared  to  those  of  hydrodynamic  motions.  A  pulse  will  then 
produce  a  hot  channel  which  is  initially  at  increased  pressure.  The 
experiments  of  Greig  et  al.  show  that,  over  time  scales  x  *  2(R1-RQ) /cg, 
the  hot  channel  will  expand  almost  adiabatically  to  reach  pressure  equilibrium 
with  the  ambient  gas.  Here  Rg  is  the  initial  channel  radius;  Rl&  R(  x)  is 
the  channel  radius  when  equilibrium  is  reached;  and  cg  is  the  ambient  speed 
of  sound.  Energy  losses  and  radiation  should  be  less  than  20%. 

Following  pressure  equilibration,  the  channels  cool  and  expand  according 
to  the  equation 


R2(t)  -  Rt2  +  4e(t-x),  (1) 

where  a  is  the  effective  thermal  diffuaivity  and  t  is  the  tine  elapsed  since 
energy  was  deposited.  The  experimental  data  show  that 


(2) 


a  -  Cg  RiA  r 

in  which  k  *  66±6.  This  gives  values  for  u  which  are  3-4  orders  of  magnitude 
larger  than  that  due  to  nonturbulent  thermal  conduction  in  ambient  air.  At 
the  same  time,  channel  turbulence  grows  to  a  fully  developed  state.  These 
observations  suggest  that  cooling  occurs  through  the  entrainment  of  cooler 
ambient  gas  by  the  turbulent  flow  field. 

Recent  theoretical  analysis  and  numerical  simulations  have  provided  both 
a  viable  mechanism  for  the  generation  of  turbulence  in  hot  channels  and  a 
useful  framework  for  analyzing  and  predicting  the  properties  of  experimental 
channels.  In  the  following  sections,  we  review  recent  progress  in  the  theory 
of  turbulent  channel  cooling.  We  show  that  any  deviations  from  cylindrical 
symmetry  when  a  pulse  deposits  energy  in  the  surrounding  gas  will  produce 
considerable  large  scale  vorticity.  lhis  vorticity  persists  and  smaller  scale 
structure  develops  as  ambient  gas  mixes  with  the  channel,  we  identify  basic 
asymmetry  classes  and  state  equations  for  the  vortex  strength  and  mixing  time 
scale.  From  the  vortex  strength,  we  then  obtain  the  effective  thermal 
diffusivity  for  comparison  to  experimental  data,  the  final  section  presents 
data  from  numerical  simulations  of  the  various  types  of  asymmetry. 

Theory  of  Convective  Cooling 

As  a  hot  channel  relaxes  to  pressure  equilibrium,  any  deviations  from 
cylindrical  symmetry  will  lead  to  asymmetries  between  the  local  gradients  of 
pressure  and  density.  This  will  cause  vorticity  to  be  generated  according  to 
the  following  equation! 
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(3) 


dC 

*)/p2  , 

where  £  =  7  x  v  is  the  vorticity,  v  is  the  fluid  velocity,  p  is  the  density, 
and  P  is  the  pressure.  All  variables  are  functions  of  the  position  r  and  the 
tine  t.  The  three  generic  types  of  asynaetry  are  as  follows: 

(1)  noncollinearity  of  successive  pulses,  (2)  pulses  with  concentric, 
noncircular  <e.g. ,  elliptical)  energy  contours,  and  (3)  curvature  or  kinking 
of  the  hot  channel.  In  case  { 1 ) ,  the  shock  wave  produced  by  the  second  pulse 
propagates  through  the  channel  produced  by  the  first.  Because  the  pulses  are 
not  aligned,  the  source  tern  in  Bq.  (3)  is  nonzero.  An  important  extension  of 
this  case  is  the  situation  in  which  the  interior  of  a  pulse  or  discharge  is 
nonuniform;  for  example,  the  pulse  or  discharge  may  contain  "hot  spots"  or 
might  have  a  filamentary  structure.  Several  subregions  will  contain 
significant  vorticity,  and  mixing  will  be  increased  over  case  (2),  in  which 
the  contours  by  definition  obey  a  similarity  transformation. 

We  have  found  that  the  residual  vorticity  for  the  three  generic  asymmetry 
classes  may  be  represented  in  terms  of  one  or  more  vortex  pairs  [2,4]  of 
strength  or  circulation  where  a  «  r,  9,  or  x.  If  we  align  the  s-axis 

with  the  average  position  of  the  beam  path,  we  find  that  <z  is  the  only 
nonzero  cylindrical  component  for  the  two-dimensional  asymmetry  c lessee  (1  and 
2)  and  that  <z  and  Kg  are  important  for  curved  or  kinked  channels  in  class 
(3). 


To  discuss  the  transport  of  energy  by  the  residual  rotational  flow  field, 
we  consider  one  vortex  pair,  aligning  the  y-axis  with  the  centers  of  the 
vortices  and  placing  the  x-axis  midway  between  them.  The  fluid  velocity  along 
the  x-axis  is 


v  (x,0)  -  k  6/ir(x2+S2)  ,  (4) 

X  z 

where  26  is  the  displacement  of  the  vortices.  Me  have  found  that  6  is  of  the 
same  order  as  the  length  scale  (radius)  Sq  of  the  local  density  gradient 
(channel).  Me  define  a  mixing  time  scale  T^.  as  the  time  required  for 
a  fluid  element  to  move  from  (x,y)  «  (-SQ,0)  to  (x,y)  -  (SQ,0). 

Integrating  Sg.  (4),  we  have 

Tmix  *  21*0  (S0/3  +  *2,/l  **1  6  *  (5) 

Me  have  computed  the  vortex  strength  for  the  three  generic  asymmetry 
classes  and  for  class  (3).  All  have  the  same  basic  tom 

K  -  (0*  T/2)  In  (p„  /pQ)  f  ,  (6) 

where  is  a  characteristic  velocity  of  the  expanding  channel  boundary 
(U  ~  c  ) ,  p  is  the  ambient  density,  p.  is  the  minimum  density  in  the 
channel,  and  t  has  been  defined  as  2(Rt  -  RQ)/cs.  The  quantity  f  is  a  form 
factor,  which  is  <  i,  the  nominal  range  being  0.15-0.33. 


\ 


We  may  now  derive  an  approximate  equation  relating  the  effective  thermal 
diffusivity  a  and  the  vortex  strength  <•  We  consider  a  pair  of  oppositely 
oriented  vortex  tubes  which  are  contained  in  a  cylindrical  region  of  radius 
Sq,  so  that  the  separation  of  the  vortex  axes  is  26  ~  2SQ.  The  z-axis  is 
parallel  to  the  axes  of  the  vortices  and  the  (y,z)  plane  contains  the  axes. 
Using  Bj.  (4),  we  find  that  the  mass  flow  of  ambient  gas  into  a  length  t  of 
the  boundary  cylinder  (and  between  the  vortices)  is 

L/2  5/2 

Pm  f  dz  f  dy  v  ~  P.  L  |  k|/»  (7) 

-L/2  -5/2 


We  integrate  over  only  half  the  separation  of  the  vortices  to  account  for 
their  finite  radii  ("■  i/2,  as  indicated  by  the  numerical  simulations ) .  From 
Eq.  ( 1 ) ,  we  compute  the  mass  flew  rate  due  to  this  entrainment  to  be 


dm 

dt 


R2(t)  -  Ri2 

P-  L  -  t  -"t - -  4ap-L  • 


(8) 


From  Bqs.  (7)  and  (8),  we  then  have 


a  ~  )  k|/4* 


(9) 


Substituting  for  <  end  t,  we  have 


«  -  cg  (Ri-R0)  ln(p./P0)|f|/4*  . 


(10) 


i 


For  an  electric  discharge  depositing  300  -  600  J/m  in  air,  we  use 


adiabatic  expansion  to  estimate  Rj  *  1.4  cm  *  3  Rg,  and  Pj/Pg  **■  10.  Comparing 
Eg.  (10)  with  the  experimental  data.  Eg.  (2),  gives  an  empirical  form  factor 
|f|  *  0.12,  which  lies  at  the  low  end  of  the  range  which  we  would  expect. 
Similarly  for  a  C02  laser  depositing  approximately  9  J/m  in  nitrogen  doped 
with  SFft,  we  estimate  R^~  1.4  RQ  ~  0.5  cm,  P^/Pg  ~  2.  Equations  (2)  and  (10) 
predict  a  value  of  |fj  ■  0.96.  this  high  value  occurs  because  the  multimode 
structure  of  the  laser  pulse  produces  several  hot  spots,  each  pair  of  which 
interact  to  produce  vortex  filaments.  To  treat  the  problem  properly,  we  would 
have  to  propose  a  distribution  of  hot  spots  which  satisfies  the  constraints  of 
total  energy  deposition  and  marginal  detectability  by  laser  burn  patterns.  We 
could  then  analyze  all  pairwise  interactions  and  compute  a  vortex  strength  for 
the  channel.  For  pulses  containing  a  few  (N  -  5)  hot  spots,  we  might  expect 
an  enhancement  in  diffusivity  (or  effective  total  form  factor)  of  'S  or  so, 
depending  on  the  average  separations,  overpresssures,  and  spatial  scales  of 
the  hot  spots. 

Numerical  Simulations 


Using  the  FAST2D  and  FAST3D  computer  codes  [5-7]  we  have  performed 
numerical  simulations  of  vorticity  generation  in  the  generic  asymmetry  classes 
and  the  case  of  a  laser  pulse  with  hot  spots.  Descriptions  of  algorithms, 
grids,  and  initial  conditions  for  the  simulations  appear  elsewhere  [4,7-8]. 

In  all  cases,  energy  deposition  was  instantaneous;  the  ambient  gas  was  air; 
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and  the  initial  conditions  were  quite  similar  to  those  in  the  corresponding 
experiments  of  Greig  et  al.  [1-3].  The  figures  belcw  show  only  the  fine  grid 
surrounding  the  hot  channel,  and  time  t  »  0  corresponds  to  the  instant  at 
which  energy  is  first  deposited.  Due  to  limited  resolution,  we  model  only  the 
large  scale  vorticity  and  cannot  follow  a  cascade  to  smaller  scales. 

Figure  1  illustrates  the  case  of  sequential  noncoLLinear  pulses  [4] .  The 
first  pulse  deposits  energy  in  the  center  at  time  t  *  0.00  ms,  while  the 
second  deposits  energy  to  the  right  of  center  with  a  displacement  of  one  pulse 
radius  at  t  -  1.00  ms.  The  peak  overpresssures  at  the  instants  of  energy 
deposition  are  4.7  and  2.4  atm,  respectively.  Figures  1 (a)  and  (b)  show 
density  contours  and  velocity  vectors,  respectively,  at  1.26  ms.  The 
noncollinearity  of  the  pulses  has  produced  a  hot,  oblong  channel,  and  after 
the  shock  wave  from  the  second  pulse  has  propagated  away,  a  vortex  filament 
pair  remains.  The  separation  is  2$  “1.1  cm  ~  R1  and  the  enhanced  flow 
between  the  vortices  covers  a  vertical  distance  of  approximately  6,  as  assumed 
in  the  derivation  of  Eq.  (9).  The  vortex  strength  falls  within  20*  of  the 
value  given  by  Bq.  (6). 

To  illustrate  asymmetry  class  (2),  we  study  a  pulse  with  elliptical 
contours  [4].  Symmetry  considerations  permit  us  to  simulate  only  the  upper 
half  of  the  channel.  Figure  2  shows  density  contours  at  t  *  (a)  36  in  and 
(b)  7.8  ms.  In  Fig.  2(a),  ws  see  the  shock  wave  propagating  away  from  a  hot 
elliptical  channel  having  an  equilibrium  area  of  ~  0.40  cm2.  Hot  visible  are 
two  oppositely  directed  vortex  filament  pairs,  one  on  either  aids  of  the  major 
axis.  Fig.  2(b)  shows  that  the  flew  between  each  pair  and  the  slow. 


horizontal,  outward  movement  of  the  pairs  have  pulled  fluid  into  the  channel 
along  the  major  axis  and  out  of  the  channel  along  the  minor  axis.  The  value 
of  a  coming  from  the  calculation  is  a  factor  of  more  than  6  too  low  when 
compared  to  experiment  [3],  primarily  because  the  interior  structure  of 
laboratory  pulses  was  not  smooth  [8]. 

Figure  3  shows  vorticity  contours  computed  by  FAST3D  at  t  *  i92ys  for  a 
ha If -wave length  section  of  a  sinusoidally  curved  electric  discharge,  which 
oscillates  (spatially)  about  the  z-axis  [71.  Figure  3(a)  shows  the  vorticity 
in  a  giver,  x-y  plane  (corresponding  to  Zz) ,  and  Fig.  3(b)  shows  the  x-z 
symmetry  plane  (corresponding  to  Cy) •  The  magnitude  of  a  agrees  well  with 
the  experimental  value. 

In  the  case  of  a  laser  pulse  with  a  circular  envelope  and  smooth  energy 
contours,  we  would  expect  little  vorticity  to  be  produced  by  the  mechanism 
proposed  in  this  paper.  In  fact,  experimental  data  on  laser  pulses  with 
roughly  circular  envelopes  have  given  a  value  of  a  -  250  cm2/s  [3].  we  have, 
therefore,  investigated  the  possibility  that  the  pulses  contain 
nonuniformities  or  hot  spots  due  to  the  existence  of  nonaxial  modes  [8].  In 
Fig.  4,  we  show  (a)  pressure  contours  at  t-0 . 00  ms  and  density  contours  at 
t*(b)  0.06  ms,  (c)  0.41  ms  and  (d)  1.24  ms  from  a  numerical  simulation  in 
which  the  laser  pulse  is  assumed  to  contain  seven  hot  spots.  Each  hot  spot 
produces  a  shock  wave  which  sweeps  through  the  shocks  and  the  density  minima 
resulting  from  the  other  hot  spots.  These  interactions  produce  vortex  pairs 
which  reside  between  the  six  outer  density  minima  which  are  left  after  the 
shocks  have  propagated  away,  ht  later  times,  the  outward  movement  of  the 
vortex  pairs  and  tha  enhanced  flows  between  each  pair  distort  the  channel 
boundaries  and  eventually  pull  the  channel  apart. 


Our  numerical  simulation  yields  a  ~  100  cm2/s,  which  is  within  a  factor 
of  2.5  of  the  experimental  value.  More  realistic  initial  conditions  should 
produce  even  better  agreement.  In  addition,  Greig  et  al.  have  studied  a 
roughly  circular  laser  pulse  containing  only  axial  modes  and  have  found  a 
significantly  reduced  diffusivity  (a)  of  3  cm2/s,  as  the  theory  would 
predict. 


Conclusions 


We  have  reviewed  a  general  theory  of  convective  cooling  of  hot  channels 
produced  by  a  beam  propagating  through  a  gas.  The  theory  provides  good 
agreement  with  experimental  data  on  laser  and  electric  discharge  channels;  the 
range  of  energy  deposition  covered  by  the  data  is  roughly  two  orders  of 
magnitude.  We  predict  that  any  deviation  from  cylindrical  symmetry 
by  beam  pulses,  including  noncollinearity,  will  produce  long-lived,  local 
vorticity.  Pulses  having  interior  fluctuations  in  energy  will  exhibit 
enhanced  cooling  over  those  with  smoothly  distorted  energy  contours.  A  goal 
of  future  work  is  to  develop  better  models  of  energy  deposition  and  the 
turbulent  cascade  to  permit  more  realistic  simulation  of  the  entire  process  of 
channel  production  and  cooling. 
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1.26  ms 

(a) 


1.26  ms 

(b) 


Figure  I 

Vortex  pair  generation  by  sequential,  noncollinear  pulses  of  energy.  T 
diagrams  show  (a)  density  contours  and  (b)  velocity  vectors  after  the  shock 
wave  from  the  second  pulse  has  propagated  far  away  from  the  channel. 


Contours  of  vorticity  produced  by  a  numerically  simulated  electric 
discharge  with  a  sinusoidally  curved  axis.  Ihe  diagrams  show  typical 
vorticity  distributions  in  (a)  transverse  and  (b)  longitudinal  planes. 


u.4i  ms 

(c) 


1.24  ms 

<d) 


Figure  4 

Pw**um  contours  and  (b)-(d)  density  contours  shoving  energy 
deposition  by  a  simulated  laser  pulse  vith  interior  hot  spots.  The  resulting 
vortex  filament  pairs  propagate  outvard,  pulling  the  channel  apart. 
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